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NOMENCLATURE 


a^ Ambient speed of sound 

a Radius of the circle in 5-plane 
C Chord of the Aerofoil 

D Drag force on the aerofoil 

F Field point in the far field 

g Forcing function defined by Eq. (2.4.7) 

i Coordinate along the imaginary axis and = -1 

K Vortex strength 

L Lift force on the aerofoil 

Mach number 

n Exponent = 3 in Eq. (3.3.16) 

p pressure in the flow field 

•4 

q perturbed velocity vector 

R Distance from the field point to the source point in xy plane 

r radial coordinate 

s source point location 

t Time 

^ Stress tensor 

u x component of the perturbation velocity 
Up Velocity of the vortex 

v y component of the perturbation velocity 
V Velocity vector in the flow field 
W Complex potential in 5~plane 

x . x coordinate in the physical plane 

y y coordinate in the physical plane 

Z ’ Complex coordinate 

Cp Pressure coefficient-defined by Eq. (3.4.3) 
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Greek Symbols 


a = K/U a 

CO 

p 

y 

6 

i 

T1 
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<p 
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Nondimcnsional vortex strength 
Circulation imposed around the aerofoil 
Strength of trailing vortex sheet per unit length 
Delta function 

Coordinate in the complex plane 
Coordinate in the complex plane 
Angular coordinate 
Circulation around the aerofoil 
Complex coordinate = ? + iT) 

Total velocity potential for the flow 
Perturbation velocity potential 
Stream function for the flow 
Kinematic viscosity 
Density of the fluid 


Subscripts 


i,j -* Coordinate indices 

0 Initial conditions 

00 Infinity conditions 

Superscripts 

> Perturbed quantities 

~ Nondimens ional , quantities 

- Refers to complex conjugates 
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I. INTRODUCTION 


This report is concerned with the acoustic field analysis asso- 
ciated with a vortex of modified strength interacting with an airfoil. 
The general theory of aerodynamic sound has been dealt with in detail 
in Ref.l'. The basic procedure followed here to determine the sound 

field for the vortex-airfoil interaction has been that due to Light - 
2 

hill o Lighthill’ s method of determining the aerodynamic sound field 

is to firs't determine the flow field from fluid mechanical principles 

and then to analyze, independently of density in the complete flow 

field, the propagation according to the principles of acoustics. 

3 

In a second paper, Lighthill applied his general theory to con- 
sider sound generated by turbulence and in particular carried out a 
detailed theoretical examination of the sound field of a turbulent jet 
at low subsonic speeds. From the analytical treatment of the aero- 
dynamic sound problem by Lighthill, it was pointed out that solid sur- 
faces even when they are rigid, might well play an important role in 
the phenomena of sound generation. But a detailed study of the influ- 
ence of solids boundaries on the mechanism of sound generation was 

4 

presented by Curie . He extended Lighthill' s theory to study the effects 
of fixed solid boundaries present in a given flow on the sound that is 
generated aerodynamically and his results distinctly show that the far 
field sound generated is now due to (i) .the quadrupole field distributed 
in the region where all the fluid fluctuations are dominant, and (ii) to 
the surface distribution of dipoles, which represents a fluctuating force 
exerted by the solid boundaries on the flow. This second contribution is 
due to the presence of the solid boundaries and as pointed out by Curle^, 
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the intensity of sound due to the dipoles on the surface becomes increas- 
ingly important at low Mach Numbers. 

It is clear from the studies of Lighthill and Curie that conventional 
aerodynamics becomes the starting point for the study of sound field, A 
particular case of some practical interest is the noise produced by heli- 
copter blades. 

The generation of lift by a helicopter in forward flight is accompa- 
nied by the development of an aerodynamic flow field of great complexity. 
Associated with the flight of the helicopter is the characteristic noise 
called the blade slap or bang. One of the mechanisms generating this 
noise is explained to be caused by the fluctuating force on the blade due 
to fie interaction with the vortex from another blade. In this particular 
case of the blade-vortex interaction, the blade interacts with the velocity 
field of the vortex but does not cut through it. This problem of the blade- 
vortex interaction has been studied by Sears'*, among other people. Sears 
has based his analysis on the classical, incompressible, two-dimensional 
unsteady airfoil theory and further assumed that the part of the blade 
that produces the noise is many blade-chords in length and that it passes 
through the field of the vortex lying parallel to and below it. According 
to this model, the lift on the blade is determined from the gust-entry lift 
function of Karman and Sears and the fluctuating lift is what is employed 
to determine the sound generated aerodynamical ly. A fluctuating force, like 
the lift on the blade, can be represented according to acoustic terminology 
as a concentrated dipole, of appropriate strength, at the origin and thus 
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with the helicopter noi^ e , by assuming a long blade passing obliquely 
over a vortex; that is, the interaction being modeled as a two- 
dimensional airfoil ±n &n oblique gust. 

The unsteady lift on the blade due to the blade-vortex interaction 
has been calculated using linear unsteady aerodynamics in a manner 

g 

analogous to the the theory of Sears developed purely for two-dimensional 
interactions. This fluctuating lift force is then employed to calculate 
the associated sound field. 

5 7 

In the above mentioned studies of Sears and Widhall , the far field 
sound due to the vortex-airfoil interaction is calculated just by deter- 
mining the fluctuating lift force on the airfoil. The perturbed flow- 
field around the airfoil, which would constitute the quadrupoles of 
Lighthill’s theory, have been completely ignored. To what extent this is 
justifiable can only be understood after undertaking a systematic and 
thorough study of all the individual contributions to the far field sound. 
Secondly in the model employed by Sears to study the vortex-airfoil 
interaction, it is assumed that the vortex is held stationary while the 
airfoil is in uniform rectilinear motion. This is a restriction on the 
motion of the vortex. A vortex free to move while interacting with the 
airfoil seems to be a more appropriate model to postulate for the vortex- 
airfoil interaction. Again, the importance of the free vortex-airfoil 
interaction can be established after a detailed analysis of the problem 
from the fluid mechanical principles. 

Therefore, with a view to clearly understand the various mechanisms 
that enter into the determination of the sound field due to the vortex- . 
blade interaction, a free vortex interacting with an airfoil in uniform 
motion is modeled and analyzed. In Ref.l classical, incompressible, 
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inviscid and two-dimensional unsteady flow has been employed for the 
analysis of the vortex-blade interaction. In that report is presented 
a complete and detailed analysis of (i) the motion of the free vortex, 
(ii) the forces acting on the airfoil, (iii) the extent of the perturbed 
flow field around the airfoil. The determination of the sound field due 
to all the acoustic sources, not only on the airfoil surface (dipoles) 
but also due to the ones distributed in the perturbed flow field (quadru- 
poles) due to the vortex-airfoil interaction, has been presented here. 

It is important to point out that since inviscid flow assumption has 
been employed for the study of the vortex-blade interaction, the quadru- 
poles that have been mentioned above in the perturbed flow field are 
entirely due to unsteady flow field. 

The effects of thickness of the airfoil on the sound radiation are 
studied by employing a symmetric Joukowski airfoil for the vortex-airfoil 
interaction. 


4 



II. ACOUSTIC FIELD ANALYSIS 


2.1 Sound radiation in a Plane 


To determine the density fluctuation at any point r , at any 
time "t" due to a fluctuating fluid flow within volume bounded 
by the airfoil surface S and the outer boundary 2 (Sketch No. 2) 
one requires a detailed knowledge of the pressure and velocity fields 
both on the airfoil surface and the surrounding region . Such 
information has been obtained in Ref.l for a two-dimensional, unsteady, 
motion of vortex Interacting with a flat plate airfoil in an incompres- 
sible, inviscid potential flow. Hence before the details of the flow 
field are used to calculate the density fluctuation p' at the field 
point F external to , Eq. (2.4.16) in Ref.l, needs to be simplified 


to the two-dimensional case. From the assumption that in region V 


0 * 


the fluid is incompressible, Eq. (2.4.16) in Ref.l simplifies to (See 
Appendix A) 

-» f 1 /- tT] 

p' = < — — - div div / dV 

l kn r 


0 


V 0 
n • [pi] 


1 /* n * 

div / 

tar a Q J s 

Pp f 

hnr a^ / 


dS 


n - U [Sv/Sx] 


dS 


kn 


0 

-0 t n • [D „ ql .. 

dxv / dS 


! 


(4.1.1) 
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The last two integrals refer to the contribution from the dipoles 
distributed on the surface of any arbitrary aerofoil held stationary in a 
medium where the fluid is in uniform motion with a velocity U . But 

CO 

specifically, for a flat-plate aerofoil they are identically zero 
because of the surface boundary condition, namely, n • V = 0 on S 
Therefore, for the flat-plate aerofoil in an incompressible flow, 

Eq. (2*l.l) further reduces to. 


>'(? ,t) = 


div div f dv 

tar X C 


kix 


div 


/ 

S ■ 


n • [p'l] 


dS 


} 


where any function [G] is defined as 


[G] = G 


H 


x + r 


O^ 1 “ *0' 


and 


( 2 . 1 . 2 ) 


- • *S )2 + (1 ■ { (y F • y s )2 + - Z s )2 } 


From Eq. (2.1.2), we observe that the far-field density perturbation 

p'(r^,,t) is due to both the volume distribution of quadrupoles and 

1). 

surface dipoles, similar to that expressed by Curie in his paper 
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except that the present results indicate the effect' of acoustic- dipoles 
and quadrupoles present in a uniformly moving medium. 

In Eq. (2.1.2), the integrands are all three dimensional quantities. 
They are functions of all the three coordinates (x,y,z). But the fluid 
mechanics problem of the vortex-aerofoil interaction, from where the 
above integrands have to be expressed, is solved as a two-dimensional 
problem. Therefore Eq. (2.1.2) has to be integrated with respect to 
Z and expressed as a function of x and y only. This is done below. 

The field point F is chosen to be in the — 0 plane (see 
sketch below) . 
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All the flow quantities like the perturbed velocity q .and the- pressure 
p are all only varying in the xy plane and are all independent of 
z . Hence the density fluctuation p' at F which is now a func- 
tion of x,y and t , is obtained by integrating Eq. (2.1.2) with 


respect to z c , from -» to +co . Consider the integral 

. O 


where 


Denoting 


-fas 

/?-.-/// 


0 


dz s j d* s dy g 


R q 


r = 


+ ( x * *0 {(y F - y s )2 + 4 } 


(Xp - x s ) 2 + (1 - i^)(y F - y s ) 2 - 


R 


it yields for 


- \/ r2 + a - >£) 


(2.1.3) 


■ Therefore the integral now becomes , 



[T] dV Q - 
' r 



[T] dr 


\/ 1 - M 2 V r 2 - R 2 

CO 


dx s dy s 


(2.1 .h-) 


In the above equation, the integral in the parenthesis, on the right 
hand side, is defined as 
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(X> 



[T] dr 
y/7 - R 2 



T(R,t) 
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r + H x \ 

M 

a 0 (i - *J/ 



R 


dr dr 



T(R^t) dr 



Note that the above integral exists only for 

R + M x ) 

because the signal that the observer at "F" receives at time 
is only that originated at 

R + M x 

CO 

a 0 (1 - 4 ) 

and for all other t , T(R,t) = 0 

Hence the integral in Eq. (2.1.5) becomes 




g- [T(R,t)] log {-2R/a (l-l^l 

' 0 1 




(2.1.5) 


»T£ff 


( 2 . 1 . 6 ) 
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where 


[T(R,t)] = T R,t = t - 


M x + R 


a 6 (1 - *0- 


Hence from Eqs . (2.1.6) and (2.1.7) , we have 


/ 


IT] dV 


0 


0 


* 0 ( l - ^) 5/2 


// 


[T(x s ,y s )] log \ dx g 


0 


where 


R 


" x s )2 + - ^Hy r - y s )S 


and 


*4. = 2R/a 0 (l.-^) 

Similarly in Eq. (4.1.2), the surface integral becomes 


/ 


n • [ P I] 


dS = 


V 1 ■- 


/-■ 


[pi] log dxg 


Therefore, the density fluctuation p '(x^,y t) is given by 


p'(x F ,y F ,t) = 


** a o (l ‘ 


' fj 

u 


div div / / [T] log R^ dx g dy, 


+ div 


A 

ivi n • 


[pi] log R x dXg 
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(2.1.7) 


dy s 

( 2 . 1 . 8 ) 


(2.1.9) 
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where 



Furthermore since the fluid fluctuations responsible for the aero- 
dynamic sound generation are solved assuming the flow in region R^ to 
be incompressible, it would seem reasonable to identify the sound 
propagation region, that is the region external to R Q 
in the xy -plane, also as one of low Mach number. Therefore the retarded 
quantities now become 

/ R 

[G] = G ft - ~ 

' \ a o 

where R ~ ^(xp x^) • 4 - “ Yg) ’ (2.1.11) 

In the region external to Rq in the xy plane, (Sketch No. 1 ), sound 
is propagating in a uniformly moving medium. In this region the pressure 

and density are related by 

p'(VV^ = a 0 (2.1.12) 
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Thus the pressure fluctuation at point F is given by, 

p'Cx^y^t) = j div div j J [T] log R x dx s dy g 


0 


K 


0 


2n a 


0 


/" 


div / n • [pi] log R x dx g 


(2.1.13) 


where assuming that M OT is small, we obtain 

[T] = T(t - R/a Q ) 


( 2.1.1k ) 


r = y (xf ” x ^ 2 + (yp " ys) " 

and ^]_ = ^ R / a o * 

2 .2 Far Field Simplification 

Any quantity [G] is defined to be a quantity G referred to 
its value at the retarded time (t - R/a^ ) • • Thus 

a 



1 

l Sr. 

•log R 1 

"dG ‘ 

Sr 

[G] log R 1 


[G] 

X 

— 


X 

' 1 

h hx i 

a o 

. dt _ 

Sx. 

X 


&X. 


Hence for low Mach number flows. 


( 2 . 2 . 1 ) 


a I 


1 1 

[G] log R n 

'dG ' 


- — • 

[G] log R 

\ r+j 


— 


Sx. j 

X ’ 

X 

1 _ 1 

* a o 

. St . 



Sr 


Sx. 

X 


( 2 . 2 . 2 ) 


Out of the two terms in the curly brackets in the above equation, for 
large distances away from the origin, that is, for large R , the 
first term becomes smaller as compared to the second term. In fact, 
this is what defines a "radiation field". Thus in the far field, 





| log R, 

'3g ‘ 

[G] log Rj^ 

1 ~ 

1 a o 

. dt 


Sr 


Sx. 

x 


( 2 . 2 . 3 ) 
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where 


dR 

dx. 


R - R 
Fi Si 


and 

Therefore we obtain 


R = 


x «) s + (y F - y s ) 2 



a 

IG] log R- 

( ( R Fi - R Si ) log R 

'ae ' 


— "f* 

— 

ax. 

1 

X. 

a o E 

.St . 


Similarly we have. 


dx. bx.. ' 
i J 


[G ] log ^ = 


( R Fi - R Sj)( R Fj - R Sj) X °S *1 

2 t >2 • 

a o R 


&G\. 

11 


at" 


2 ^i - V^l - V 
— o R 


ac. ; 

hi 


L at 


Thus Eq. (2.1.15) can now be written as 


p'tvv^ 



(R - R )(R - R ) 

^ Ti Si A Fi Si ; 

a 2 R 2 

a o R 


2(R - R )(R - R ) 

k Ti Si A Fi Si ’ 


r a 2 T. . 
11 

L at 2 j 


log R 1 


a o R 
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/ n i 


(R - R ) 

v Fi Si ; 

a o R 
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at 
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(2.2.4) 

(2.2.5) 


( 2 . 2 . 6 ) 

dy S 

(2.2.7) 



Note that all the quantities in the above equation are referred 
to the "retarded time". The pressure perturbations measured at the 
field point F (coordinates ) at the instant of time 't' are 

actually the ones emitted by the source at S (coordinates ) at 

an earlier time 

\\ - «*! 

t 

a 

0 

where a is the uniform speed of sound propagation. The definition 
of the term "far field" always refers to points , whose distance from 
the source, - itgl is much greater than the largest wavelength of 

the propagating sound wave. 


F 



SKETCH No. 2 
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Referring to the sketch above., let refer to some maximum 

• - * 
characteristic dimension of the region . Consider a point in 

the far field F where the observer is located and the distance 

, jr^l is much greater than . Consider two source points at 

S 1 and Sg within region . Retarded time for the source 


is 



Hence 


At 





Therefore 



( 2 . 2 . 8 ) 


This means that the time difference in the two signals reaching the 
observer is at best equal to the time taken for the signals to travel 
the' maximum characteristic distance in the region 
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Consider now the time scale of the fluid fluctuations in the 
flow field in . The flow velocity is V . This can be approxi- 
mated to be the free stream velocity, for discussion here. Hence 

the time scale of fluid fluctuations is of the order L,/U 

CO 


Time scale of fluid fluctuations 

in V. = L-./U 

0 V ® 


( 2 . 2 . 9 ) 


Thus 


— „ 


( 2 . 2 . 10 ) 


Therefore, if the Mach number of the flow is low (M^ « l) , all the 
fluctuations in region will more or less reach the observer at 
the same time. Thus for low Mach number flows it seems reasonable to 

ignore the retarded time. Hence Eq. (2.2.7) now becomes 


p'OvVt) 




*1 


2 J2. 

a o E 
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Defining the following non-dimensional quantities 


P' = 


i p, 


. T = Hr- ; R, = 2 M R ; 

•13 i„ B 2 1 


R 


R = a , 


0 “» - ■ 0 CO 

U t 

^ CO 

t = 


z P c 

and a - C/4 




and expressing p / (x FJ y Fi ,t) = Py + Pg 

we have , after substituting in Eq. (2.2.11), the following; 

• 2 - ^ \ . c? /^2 k 


(2.2.12a) 


, £ ( h %i) 

? v (VV fc) 2n JJ \ R 2 \ at 2 / 


r*-F 


2 (K - R )(&. - S- .) /2ff. A) 

^ ^ . 1 ^jjdx s dy s 


R 


(2.2.12b) 


and 


M 2 f (Rp. - ll.) log /^P \ - .. . 

<vV> - t t)" 3 ' 


From Eq. (2.2.12), the term p' r refers to the pressure propagated 
into the far field due to the unsteady fluid fluctuations T. . within 

xj 

2 

the region . Then this is what Lighthill identifies as the 

quadrupole effect. Similarly p' is due to the distribution of 

u 

h 

the singularities on the surface. As Curie points out, these are 
surface distribution of dipoles and thus Pg is then the dipole 
effect. 
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Furthermore consider the term p' in Eq. (2.2.12). It is 

D 


M 2 ' 


^( X F ,y F ,t ) “ 7 /" i 

27T j 


a?\ 

r 108 R i W dXs : 


R . - R . 
Fi Si 


Considering the maximum characteristic dimension of the aerofoil 

and comparing it with the location of F , the field point, it is 

in order to assume R,,. «R_.. . Therefore, now 

Si Ti J 




K; ^i lo * E Pi r (%\ 


where 


27 T 


R F- 


Hw) 


for- 


< Si g, ~ 

#v v ro 

2tt Rp at . 


S(S - / 


P dx. 


(2.2.13) 


(2.2.14) 


= total resultant force exerted upon the 
fluid by the solid boundaries 


*fi 


= 2M 


h 


> *F “ *F + y F 


(2.2.15) 


and n^ the surface normal is independent of time. 


Therefore 3P /St is the time rate of change of the resultant force, 
which from acoustics represents an equivalent dipole . As mentioned 
in the earlier chapters, the study of sound field just by the knowledge 
of the time varying force alone, completely ignores the effect of the 
quadrupoles in the region R Q . The validity or otherwise of such 
an assumption can be established by computing the sound field from 
Eq. (2.2.12). 
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4.3 Computation of the Sound Field 


(a) In Ref .1, where the analysis of the flow field is described, the 
forces acting on the airfoil, namely, the lift and drag forces, due 
to the motion of the vortex, have been calculated and shown in Figs. 1 
and 2. As is pointed out in the preceding section, from an acoustic 
standpoint, the time varying forces can be easily represented as dipoles 
of appropriately varying strength. Hence the pressure fluctuation pro- 
pagated to the field point F due to these fluctuating forces, is given 
from Eq. (2.2.13) as 


^ M 2 S log Rp dp. 

- ± — ~ 1 ^ 

2tt Rj 5t 

where P^(t) = time varying force. Hence the pressure signal, due to 
the fluctuating lift and drag forces, is 


p'('^,y F ,t) = 


M 2 

CO 

log R Fl 

L 35 

dL 

27T 

ru 

V 




( 2 . 3 . 1 ) 


From Fig. 3, it is easily observed that 


dD dL 
— << — 
dt dt 

Therefore neglecting dD/dt with respect to 
the following simpler form , 


dL 

— , we obtain from Eq. (2.3.I) 
dt 


p's<vvt) ~ ^ /s: 


27 r 


•V 


dt 


Expressing in polar coordinates. 


„ ^ M 2 • /dL 

~ sin 0 F lo 8 ^ f • 


27T 


'l\dt 


( 2 . 3 . 2 ) 


(2.3.3) 
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The lift force L has been computed for the non-dimensional vortex 

strength, a = 1.0, for a flow velocity of = 200 ft/sec past 

the flat-plate aerofoil of chord C s 2 ft. At sea level conditions, 

this corresponds to a Mach number flow of M ro = .18. The field point 

F is chosen 25 chord lengths away so that = 100. From Eq. ( 2 . 3 . 5 ) 

it can be seen that p' is a maximum at 0 = 7r/2 . The time variation 

s c J 

of pg for the above situation is shown in Fig. 4 . It is observed 

that a peak value of this p' at 0 p = tt/ 2 occurs at t = 50 , 

that is, at a time when the vortex has moved right above the aerofoil. 

The important result to be noticed is that the free vortex motion predicts 
nearly twice as much peak value for p' as the constrained vortex 
thereby, emphasizing the importance of the vortex-aerofoil mutual 
interactions . 

(b) From the details of the flow field computation shown in Fig. 5 , 
for the given value of a = 1.0, .it is found that all the fluid fluctu- 
ations, have decayed to the so called acoustic level at a distance of 
approximately 15 chord lengths away from the origin. In fact the value 

of Cp, at 6-p — 7 t /2 and R^, = 60 (15 C) is about 0.05 • Therefore for 

sound field computational purposes, the region R^ is taken to be a 
circle of radius equal to 15 chord lengths (R^ - 60 ). Thus, shown in 
Fig. 6 , is a detail of the flow configuration adopted for computational 
purposes . 

One important thing to mention here is, that our analysis, the 
vortex considered, is a point vortex. It is well known that a point 
vortex gives rise to a singularity at its center. Therefore, in our 
problem, the vortex has to be isolated .from the rest of ■the. flow .field 
by a small circle of radius 


- 20 - 



The effect of the vortex on the flow field is felt strongly when 

it is closest to the aerofoil. From the trajectory analysis, it is 

observed athat at t e JO , the vortex is at its closest location to 

the aerofoil. The absolute value of the fluid velocity V , around 

the neighborhood of the vortex, has been computed from Eq. (3. 3*28) 

for various values of and shown in Fig. 7 . From these results, 

~ e 0 

a choice of €q — = 1.0 was made for computational purposes?, 

a 

(c) From Eq. (2.2. 12), the pressure contributions p' and p' are 

V D 

evaluated numerically in the compex (; plane as follows. The integrands 
are expressed in terms of the flow field details in Appendix D. At 
time t = 0, the vortex of strength a = 1.0, is at f = - 50 and 
■*1 q ** + 2. All the flow quantities like the velocity and pressure 
are evaluated, as explained in the earlier chapter. These quantities 
are then employed to evaluate the double integral p' and the line 
integral p' . First of all the field point F is taken to be at 
— 0 and Rp = 100. The double integral p^ is first evaluated 
for a given R g and varying from 0° to 360°, the value of R g 

ranging from R =1.0 to R = 60. Then the value of R is changed 
and integration performed with respect to Q . The quadrature sub- 

O 

routine is employed to perform the integration with respect to 0 

b 

whereas the "trapezoidal rule" is used to evaluate the integral with 

respect to R . Similarly p' is evaluated for R = 1.0. There- 

o s s 

fore we now have p/ r and p' evaluated at t = 0 , for 9„ = 0 and 

V s F 

m <v 


Rp = 100. At time t = 1.0, the vortex advances to a new location and 
again the above integrals are all evaluated and the procedure repeated 
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until t = 100 is reached. Now for the same R^ = 100; but with 6 ^ 

o ^ 

changed to 0^ = JO , the integrals are evaluated for t = 0 to 

t = 100. Thus a detailed computation of the time variation of the 

pressure fluctuation at different angular locations of the field point 

F , at a radius R^, = 100 is performed and presented in Figs. 8 (a) 

and 8 (b ) . 

(d) The pressure signal due to both the quadrupole distributions in 
the region R.^ and the surface distribution of dipoles on the aero- 
foil j which has been computed in the previous section, should be expressed 
in the more faimiliar terras in acoustics, namely the intensity of sound 
X at the field point F and its directional pattern for varying 
locations of F „ From a physical point of view, energy is being 
continuously propagated by a traveling wave and therefore intensity I 
of a traveling wave is then the "time average rate at which energy is 
transported by the wave per unit area," across a surface normal to the 
direction of propagation. More briefly, intensity is the average power 
transported per unit area normal to the direction of propagation. 

To recapitulate, it is our aim to compute the intensity of the 
sound field due to the vortex-aerofoil interaction in the space fixed 
coordinate system (Sketch No. l) in Ref.l where the airfoil is in motion 

in a fluid at rest. From acoustics it is known that the intensity vector 
I and its relation to the fluid flow quantities is given by 

* (p /S > „ 

iCRp) = n (2 .3 A) 

p o a o 
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where 


1 


and 




£t 




TT f ^ 


,t) dt 


pressure fluctuation in the flow 


constant speed of sound 


constant density of the fluid 


unit vector normal to tjde wave front. 


( 2 . 5 * 5 ) 


(2.3.6) 

I (Rp) « <p' 2 (Rp)> • (2 .3.7) 

rJ 

The pressure signal p' in the far-field^ which is 
p' s= p^ + pg , has been computed earlier for different locations 
of the field point F , for t = 0 to t = 100. These results are used 

/v» 

to determine the non-dimensional intensity I 

rs/ 

Intensity I in the far field for = 100 and different angular 
locations 0^, is calculated and shown in Fig . 10 . 

Because of the large range of intensities over which the ear is 
sensitive, a logrithmic rather than an arithmetic scale is chosen to 
define the “intensity level". Therefore the intensity level, "IL" of 


Equation ( 2 . 3 .^) can now be written as. 




I (Rg,) = ' — — I n watts/ 


2 


cm 


l^a 


0 


where 
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a sound wave is defined by the relation 


IL = 10 log L0 (I/I Q ) (2.3.8) 

where 1^ is a reference intensity. This is taken to be 10 watts/cm 
corresponding to the faintest sound wave that can be heard. Intensity 
levels are expressed in decibels abbreviated as "dB" and shown - in Fig. 10a. 

From Figs. 8 to 10, where the pressure in the far field and the 
corresponding intensities have been shown for the vortex-flat plate 
aerofoil interaction, it is seen that the maximum sound intensity in 
the far field occurs at e p = n/ 2 j that is ; ri S ht when the vortex 
has moved on top of the aerofoil. The details of the pressure signals 
at ©j, = tt/ 2 are seen in Fig. 9 . The volume distribution of quadru- 
poles contributes nearly 18$ of the pressure from that due to the sur- 
face dis.tribution of dipoles; this corresponds to a 25 $ increase in 
the associated sound intensity, as shown in Fig. 10. The sound inten- 
sity with only the concentrated dipole at the origin representing the 
total forces acting on the aerofoil agrees with the results obtained 
considering only the distribution of dipoles on the surface. These 
results can be seen from Figs. 9 and 10. The effect of the quadrupoles 
is significant in the sense that it has altered the directivity pattern 
of the sound intensity in the far-field (Fig. 10a). In terms of intensity 
levels, from Fig. 10b, the dipole distribution on the surface corresponds to 
112.3 dB vhereas the cumulative effects of both quadrupoles in R Q and the 
dipoles on the aerofoil surface correspond to II 3.25 dB at 6^, = 7 t/ 2 
for Rp — 100 and = .18. 
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Once again all the results were obtained numerically with the 


IBM 360 
for the 
to both 


computer. It took a total of 5*5 minutes of execution time 
computer to evaluate the sound intensity in the far field due 
the surface dipole' and volume quadrupole effects. 
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III. THICKNESS EFFECTS ON THE FAR FIELD SOUND 
3.1 Flow Field Study 

Up till now, the fluid mechanics of the interaction of a vortex 
with a flat-plate airfoil and the associated sound field have been 
analyzed. The thickness of the body gives rise to additional mass 
fluctuations in the fluid flow. As is familiar to all aerodynamic ists, 
such types of mass fluctuations in fluid flow are represented by simple 
sources and sinks. For closed bodies such as aerofoils, the net source 
strength representing them in the fluid flow is zero. Hence the source 
and sink distribution replacing the finite, closed body appear like di- 
poles to an observer in the far field and therefore the far field sound 
radiation is dipole-like. As before, the flow field analysis from fluid 
mechanical principles, is conducted first. 

In order to specifically illustrate the effect of thickness on the 
far field sound and with a view to utilize the complex variable technique 
that has been employed to study the interaction of the vortex with a 
flat=plate aerofoil, (Ref.l) a symmetric Joukowski airfoil is employed. 
The configuration of a vortex of strength K in a uniform stream U , 

CO 

flowing past a symmetric Joukowski airfoil of maximum thickness ratio 
e = h/C , where h = maximum thickness of the airfoil and C is the 
airfoil chord, is shown in Fig. 11. From the Joukowski transformation 
defined by Eq. (3.3.2) in Ref.l, the symmetric airfoil in the z-plane 
is transformed to a circle of radius a ; from Fig. 11 it is seen that 
the center 0^ of the circle M in the complex £-plane,is on the Q 
axis but displaced by a distance e from the origin of the 
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coordinates. The radius a^ of the circle M is slightly greater 
than the fundamental length "a" and the circle M , represents the 
symmetrical aerofoil in the complex £ -plane. Writing a^ = a(l + e), 
where e is a small quantity, we obtain for the chord of the aerofoil 
ka(l + € ) . As e is a small quantity e « 1 and thus, to a 

close approximation, the chord of the aerofoil is C = 4a and the 
trailing edge in the complex plane is represented by ^ = +a or 
^ . It is also known that the maximum thickness is at the 

quarter chord point and the maximum thickness ratio is h/c = 1 . 299 € « 

The complex potential W(i^) , for the symmetric aerofoil-vortex '' 

configuration in a uniform stream of , is given by 




D « ( £ -i + ^) + a los(c i ' W 


- iK log 


H.) 


+ iK log t >1 


+ i J l(^) log(^ - 5 X ) d^ 


a l 

a 


“ 1 /* 7(|^) l°g - J ^1 
J Q ' 


( 3 . 1 . 1 ) 


The complex velocity V(j^) is then given by 


v li'S 0 i> K > u ») “ j u „ (' - 


iK 


iK 


Ci - C01 h ' A /I 


01 


+ ^ + t r y ^i^i _ i r y ^i^ d ^i 

£1 J» ki - ?i } K <£1 - 4 >\ 


rUi ) ^ 


( 3 . 1 . 2 ) 


where = £ + ae 
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Now the Kutta condition at 'the trailing edge, £ = +a or = +a^ , , 

demands that the velocity V(£^ = +a^) be zero, at every instant of time. 
Defining the following expression for the strength of the trailing vortex sheet 

,<0 

y(5,t) = A[^ 0 (t)] 

or 


= l + ae (3.1.3) 

and satisfying the Kutta condition at 5 = +a or ^ “ +a ^ ? one 

obtains for 


7(5v fc )= A ^ .(*:)•]! 


(S 1 ~ a l } 

+ a(l “ e)p 




A[a ol (t)] « 


-ha?K 


1 - e 


+ ” 2a i 


1 < l + e > 2 (Sox - a i>^oi - a i> 


(3 .1.4) 


so that 


{ . , 2 _ ^ ^01 + ^01 " 2a l ^1 “ S l ^ 

(1 + e) (£ 01 “ a^Ctoi ~ a l^ [^1 + ‘ e ^] 

(3.1.5) 

where $ 0l = £ 0 l (t) • 

Substituting for 7 ( 5 ^, t), in Eqs . (3.1.2) and (3.1.2) one obtains for 
the nondimens ional complex potential W({^) and the complex velocity 
v(£^) the following expressions: 


^ <v rs/ , 

W ^ 01^*1 


h + 


~) + ia log (l - X t ) - ia log /f - — 

U \ ioi 


+ to log ^ - il 4 Q 1 ”~2 ■ ^ 01 (l -I&) 

1 hi*) 2 (Coi-^^or 1 ) 5 J ■ 

(3.1.6) 
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where 


I- = 


X e " 1 + 2 ti ~ 

+ log (1 + € ) 


i2(^ + e) 2(^ + e) £ 

(£ x " 1 f I?! " l| 


2(1 + e)^ + e) 


*6 = 


(?! ' 1 f log I?! " ll ?i@ ~ 1) + 2 ?i 

log (e/l + e, 

2(1 + e)(6^ + IF 2(6^ + 1) 




2e(e£. + 1) 2e 


e - 1 

+ -zg- log (1 + e) 


■?} • 


and 


e = 


1 - e 


1 + e 


a 


; and 0!^ - 


1 4 - £ 


Also 


rJ j r»# /v | 

v Ui4 0 i^ii- 


ia.. 


1 _ __ j + 


1 a i^01 


rv» /s/ 


ii ■ toi ^oi “ 1 


ia, 

+ ila 

?1 


1 - e 


^01 * £pl “ 2 

'(i,,) 2 (?01 ■ 1)( Coi ' x) 


<w 


where 


*7 " 


Cl 


~v2 log 


2(1 + t %f 2(? 1+ ?) 2 (?! ♦ ?) 


Cl 


-1 


1+e 


(3.1.7) 


(3.1.8) 


(3.1.9) 


(3.1.10) 


(3.1.11) 
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(t! ~ 1) 1 . (ti - 1) 



2(1 + + if 2?^ + l) 2 ‘ (H X + l) 3 108 

1+e 



(3.1.12) 


and as before 


- ‘ 1 - € 

e “ 1 + e 


»v /v 


The velocity of the vortex^ at , is calculated as 

before (Eq. ('3 *3 *25)) in Ref.l, and it is given by 


U 0^0l^ G l^ = 


ia. 


l - 


ia l ^01 ^ “TL 

rss T" 


Hi' loiloi-y Hi 


- ito. 


1 - e 


+ - 2 


^01 + ^01 


(1 + e) (t, 0 l " ^^oi " 


(*n “ 


I 10^ 


(3.1.13) 


where 


I« = 


&01 


- 1 


1 , ^01 ’ 1 - 

+ f log 


9 | 2 ( 1 + 6 )(t 01 + e) 2 2(t 01 + cf (? 01 + e)5 


01 


=*01 


^01 


- 1 


1 + e 


(3.1.14) 


and 


^10 = 


^01 


- 1 


^01 


>1 


2(^ 01+ l) 2 (e +1) 2e(^ 01+ l) 2 (^ 0l+ l)' 


log 




1 + e 


(3.1.15) 
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The circulation developed around the aerofoil is given by 


00 

r (t) - -J y^^d^ 


Substituting for r(§^t) from Eq. (3.1.3), it yields 

r(t) - t - (— )( ■ *>- --} 

\l - e/{(r nl - 1 )(t nl - 


(3 .1.16) 


with V = r 0 i (f) • 

With the velocity field thus being determined, the trajectory of 
the vortex and the lift force acting on the aerofoil are calculated for 
a typical case of e = 0.10 . As before the strength of the vortex is 

0C^= 1.0 and the vortex is located at t '= 0 at = (-50 + i2) 

The trajectory of the vortex and the lift force experienced by the 
aerofoil are shown- in Fig. 13 . A comparison of the lift force experi- 
enced by the symmetric "aerofoil with that acting on the flat plate 
aerofoil is shown in Fig. 14 . Also the- circulation around the aero- 
foil for the two different cases is shown in Fig. 12. 

It has been noted earlier that the thickness effect or finite- 
ness of the body can be represented by suitable distribution of sources 
in the flow. It was also observed that the net effect of such sources 
replacing the body is dipole-like at large distances. From Eq. (2.1.1) 
far field pressure is given by 


p'tv 6 ’ = is 


div div 


f dv o + div / ^ 


Ml 


dS 


'0 




n • [U (3q/3x)] 


dS - p 


/*n ♦ [?J q] 

Q div / dS 

-4 :.r * 
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Upon reducing the above equation to the two-dimensional case of sound radia 
fcion in a plane, and by integrating with respect to z (as in the flat- 
plate case) it yields for the pressure p' the following expression; 


p'C^t) = 


2*ra.(l 


^ 5 ^ 


O v co 


div div / [T] log dR^ 


,iv /" 


+ div / n • [pi] log R^ dS 


" P n f n • U 

Of CO Sx 

J 


log R^ dS 


P.0 div 


iv / n *, [tMq] log dS 


(3.1. IT) 


where, as defined earlier. 


R = 


(x F ' X S )2 + (1 “ ^ )(y F ' y S )2 5 R 1 = 2R / a O (l -*0 


and 


M^x + R \ 

a 0 d - >£)/ 

For low Mach number flows, Eq. (3.1.17) 'simplifies to 
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p'(Vt) « 


27Ta^ 


l 

R o 


div div / [T] log dR^ 


+ div 


iv /" 


fpl] log R x dS 


- p 0 / " • k i] los K i ds 

s 

“ Pq diV^ ^ ' '1^00^] i0 S dS 


* S 


(3.1.18) 


where 


R ~ 



f + (y F “ y s } ; R 1 = 2R / a 0 


‘ " and 





M + R 



) 


From the above equation, the pressure in the far field and thus the 
sound intensity can be determined, as before. 


3 .2 Far Field Sound Radiation 

In the preceding chapter for the case of the vortex interacting 
with the flat-plate aerofoil , the analysis of the acoustic field 
clearly illustrated the significant role played by the quadrupoles 
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distributed in the volume Vq enclosing the perturbed flow field. It' 
is our purpose now to study specifically the effect of thickness on the 
radiated sound to the far field. Therefore we direct our attention only 
to the surface integrals in Eq. (3.1.18). The first of the three sur- 
face integrals, as seen from the earlier chapters, is the effect of 
the dipoles distributed on the surface of the aerofoil. Even though 
this integral gives different contribution from that for the flat- 
plate aerofoil, one can specifically identify the other two surface 
integrals as definitely due to the thickness of the aerofoil. There- 
fore, grouping these three surface integrals together and writing them 
in their -nondimens ional form it yields 


Ps(V c > = 



log (2 R)ae 


J R„. “ R c . /8q.\ 

-+ 2 / — — (" / cos 0 lo S (2 M R) d0 

J Q r Vat / 


2tt 


+ 2 


- Rp - R / Cq dq\ 

f — (cos 0 — ~ + sin 9 ) log (2 H R)d0 

J r . R ' ot dt / 




(3.2.1) 


The velocity field components q , q and the pressure distribution 

x . y 

p have all been known from the flow field analysis in the earlier 
v section. As in the flat plate case, the above integral is evaluated 
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for the initial location of- the vortex £ = (-50 + i 2 ) at t* = 0 
and for all the subsequent locations of the vortex, as it moves along. 

The results for the pressure in the far field for = 100 and for 
different angular locations 0 „ of the field point are all shown in 
Figs. 15(a) and 15 (b) for the specific vortex strength a = 1 .0 . 

Once again the maximum contribution to the pressure history with time 
occurs just when the vortex is above the aerofoil. The pressure is a 
maximum for the field point location of 0 = tt/2 for R^ = 100 . 

In Fig. 16 is shown a comparison of the pressure variation with 
time, for the field point location of 0 = tt/2 and R^, = 100 com- 

puted from a single concentrated dipole model with the one computed 
from the surface distribution of dipoles. There is a slight difference 
in the results unlike for the case of the flat-plate aerofoil. This 
difference can be easily explained from Eq. (3.2.1). The first integral 
represents the concentrated dipole whereas the effect of thickness 
is eoniming from the other two integrals. 

In Fig. 17 is shown the comparison of the pressures determined for 
both the symmetric and flat-plate aerofoil cases, as predicted by the 
lift force action on the aerofoils. The effect of thickness on the 
far field is quite evident. It is clearly seen that a 10$ thick symmetric 
aerofoil has increased the peak pressure in the far field by nearly 50 $ 
over that predicted for the flat-plate aerofoil. The sound intensity, 
being proportional to the square of the pressure perturbation, is nearly 
twice as much for the 10$ thick symmetric aerofoil as is for the flat- 
plate aerofoil. This is shown in Fig. 18. The corresponding intenstity 
levels are shown in Fig. 19. 
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VI. SUMMARY AND CONCLUSIONS 


A study of the aerodynamic sound produced by a vortex passing 
past an airfoil in uniform motion is studied. 

Following theprocedure proposed by Lighthill to study how 
aerodynamic flows produce noise, the complete flow field is broadly 
divided into two separate regions, namely (i) a region of sound 
generation in the neighborhood of the airfoil, and (ii) a uniform 
madium at rest where sound is propagated due to the fluid fluctua- 
tions in the region around the airfoil. 

The vortex-airfoil interaction problem is formulated as a 
vortex of strength K being released ahead of a flat-plate of chord C 
at a specific location, at time = 0 . Classical two-dimensional 
incompressible, inviscid, unsteady potential airfoil theory has been 
employed to determine (i) the trejectory of the vortex, (ii) the 
forces acting on the airfoil, and (iii) the complete flow field 
around the airfoil. This analysis is presented in Ref.l. 

The computation of the sound field due to such a fluid mechanics 
problem is done in two ways: (i) the forces acting on the airfoil, 

specifically the lift force, from acoustics, is represented as a di- 
pole at the origin and its far-field sound radiation determined; (ii) 
not only the forces acting on the airfoil but also the fluid fluctu- 
ations in the flow around the airfoil are taken into account to com- 
pute the far-field sound. The analysis does indicate that the fluid 
fluctuations in the neighborhood of the airfoil contribute nearly 25 
to 307. to the intensity of sound in the far-field. 

The effects of thickness of the airfoil on the sound radiation 
are studied by employing a symmetric Joukowski airfoil interacting 
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with the free vortex. The increase in the sound intensity in the 
far-field is nearly twice as much for a 10% thick symmetric airfoil 
as it is for the flat plate case. 

All the numerical computations were' performed on IBM-360. It 
took less than a minute of comuter time to compute the trajectory 
of the vortex and the free forces acting on the airfoil, whereas- the 
computation of the directivity pattern of the sound intensity in the 
farfield due to both the surface dipoles and volume distribution of 
quadrupoles took a total of 5.5 minutes of execution time. 

Conclusions 

From the results of the acoustic field analysis of the vortex- 
airfoil interaction, the following conclusions can be drawn; (a) the 
estimation of the sound generated just by the knowledge of the lift 
force acting on the airfoil, underestimates the intensity in the far- 
field. The distribution of quadrupoles in the region around the air- 
foil makes significant contributions to the sound intensity in the 
far-field and thus cannot be neglected. Our earlier statements that 
the far-field sound is due to the cumulative effect of, not only the 
dipoles distributed on the surface of the airfoil, but also the qua- 
drupoles distributed in the perturbed flow region around the airfoil 
is well borne out by the results; (b) The effects of thickness on 
the sound intensity in the far field is illustrated by considering 
only the dipoles distributed on the airfoil surface. It is observed 
from the results that a 10% thick symmetric airfoil nearly doubles 
the sound intensity in the far-field over its value for :a’ flat plate 
airfoil. 
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APPENDIX A 


A.l Consider div(pq q+ p£) in Eq. (2.4.16) in Ref. 1.. From 
and momentum equations, Eqs. (2.4.4) and (2.4.5) in Ref 
equation can be written as 

i^t + ^oo * + div P? Q - Vp , so that 

div(pq q + pi) = - + U ro * V ) 

But ^ = lL + q ; U “ constant 

CO 


Hence for the incompressible flow considered. 


dlv(pq q + pi) -P 0 (^+fi„ -'v)7 


- -Kl + poi • 


Therefore n * div(pq q + pi) 


3^ 


- -Po S - (5.-^ 


But n • ^ = 0 on the body surface. 


I. 
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[div(pq q + pi)] 
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J* it - ?, 
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the continuity 
1 , the momentum 


(A. 1.1) 


(A. 1.2) 


(A.1.3) 


(A. 1.4) 



A . 2 Consider the term pV V 
Now V = U + q 

CO 

Hence p $ $ = p $(tJ + q) 

CO 

= p^ ?J + p^ q 
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3 ^ U + o 5 Q 
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00 oo 
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because n * ^ = 0 on S 
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APPENDIX .B 


B.l From Eq. (3.5.4) in Ref.l, the non-dimensional perturbation is given by 


P' = 


(p-pj/|p u !) 




= 1 - V - 2 


St 


- 1 - V V - 2 Real [ — 


/SW\ 
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(B.l.l) 


Hence the time derivative of p' is giv.en by 

. 2 - 
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= - 2 Re 
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(B .1.2) 
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where 


(B.l. 3) 
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= Uq ~ iv^ - complex vortex velocity at £ = ^ 


and 


U = Complex conjugate of U, 


Similarly we obtain for. 
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and 
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B.2 Again from earlier Eq . (B.1.2)^ we have 
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But from Eq. (b. 1.3)> we obtain 


(B.1.5) 


(B .2.1) 
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B«3 Consider the term 
>2 


br 


%c + P XX^ in Eq * ( 2 -3-6) 
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where q^ = Real part of the complex perturbation velocity 


v'(0 


qf y = Imaginary (V'(^) 


Hence we have, 
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and 



(B.3.3) 


Also' note that, 


7T — x component of — — 

a ? 2 a ? 2 


Similarly we obtain for 



fB.3^) 


(B .3.5) 
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C = ka 

Fig. 2 -Drag on the plate vs vortex location. 
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Fig. 4 --Pressure in the far field due to the lift force alone 
for the case of constrained and free vortex motion. 
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FIG. 











Fig. 12 -Circulation around the symmetric aerofoil as a function of time. 




Fig. 13 -Lift variation on the symmetric aerofoil with respect 
to t ime . 



Fig. • -A comparison of the lift forces for the vortex-flat 

plate aerofoil and the vortex-symm. aerofoil. 











Fig. 17 --Far field pressure due to a concentrated dipole — a 

comparison of the vortex-flat plate and vortex - symm. 
aerofoil interactions. 
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Fig. .19 -Intensity level (in deciBcls) due to surface distribution of 
dipoles for the flat-plate and symmetric aerof ils. 



